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Abstract. The atmospheric Cerenkov imaging technique 
has been used to search for point-like and diffuse TeV 
gamma-ray emission from the southern supernova rem- 
nant, W28, and surrounding region. The search, made 
with the CANGAROO 3.8 m telescope, encompasses a 
number of interesting features, the supernova remnant 
itself, the EGRET source BEG J1800-2338, the pulsar 
PSR J1801-23, strong 1720 MHz OH masers and molec- 
ular clouds on the north and east boundaries of the rem- 
nant. An analysis tailored to extended and off-axis point 
sources was used, and no evidence for TeV gamma-ray 
emission from any of the features described above was 
found in data taken over the 1994 and 1995 seasons. Our 
upper limit {E > 1.5 TeV) for a diffuse source of radius 
0.25° encompassing both molecular clouds was calculated 
at 6.64x10"-^^ cm~^s~^ (from 1994 data), and interpreted 
within the framework of a model predicting TeV gamma- 
rays from shocked-accelerated hadrons. Our upper limit 
suggests the need for some cutoff in the parent spectrum 
of accelerated hadrons and/or shghtly steeper parent spec- 
tra than that used here (—2.1). As to the nature of 3EG 
J1800— 2338, it possibly does not result entirely from 7r° 
decay, a conclusion also consistent with its location in re- 
lation to W28. 

Key words: Gamma rays: observations - supernova rem- 
nants: individual: W28 



1. Introduction 

Supernova remnants (SNRs) have long been thought to be 
the dominant source of cosmic rays (CR) at energies below 
100 TeV (for a review see e.g. Blandford & Eichler |1987D . 
SNR, via the diffusive shock process, are able to acceler- 
ate electrons and hadrons and meet the energetics of the 
observed cosmic rays. The TeV gamma-ray flux predicted 
from SNR is the most accessible tracer of CR acceleration 
and its detection would be convincing evidence for the 
SNR origin of galactic CR. Models of the TeV gamma- 
ray emission from SNR predict distinct spectral features, 
according to the hadronic and/or electronic nature of the 
parent CR accounting for the gamma-ray flux (see Drury 
et al. |1994 Naito and Takahara [l994| . Baring et al. |1999 , 
and references therein for a summary). 

Ground-based surveys of SNR at gamma-ray energies 
(TeV to PeV) have been carried out on several promis- 
ing northern hemisphere candidates (e.g. IC443, Tycho's 
SN R, W 51, W44, 078^2+2.1). The Whipple (Bu ckley et 
al. |1998|, Lessard et al. |1999D, HEGRA (Hess et al. |1997| at 



TeV energies, and Prosch et al. 1996 at multi-TeV ener- 
gies), CAT (Goret et al. pJ99| ) and CYGNUS (Allen et al. 
1995) groups have reported upper limits. Recently how- 
ever, the HEGRA has seen marginal evidence for TeV 
gamma-rays from the young SNR Cas-A, after deep ob- 
servation (Piihlhofer et al. 1999). In the southern hemi- 
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sphere, the CANGAROO has reported the detect ion o f 
TeV gamma-rays from SNR SN1006 (Tanimo ri et a l. |l998| ) 
and SNR RX J1713. 7-3946 (Muraishi et al. ^000|) , and if 
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confirmed, will be strong evidence in favour of the pro- 
duction of cosmic rays electrons in SNR. 

W28 (also SNR G6.4-0.1 from Green [l99^ ) is consid- 
ered an archetypal composite (mixed or M-type) super- 
nova remnant, characterised by a centrally filled X-ray and 



shell- like radio morphology (Rho & Petre 1998, Long et al. 
199ll ). TheROSATX -ray emission appears best explained 



by a thermal model (Rho et al. 1996) although Tomida et 
al. ( |1998D from the analysis of ASCA data, has suggested 
the presence of a weak a non-thermal component in the 
south west region. The limb-brightened radio emission (20, 
6 &: 2cm) shows a synchro tron s pectrum of varying spec- 
tral index (Andrews et al. |l983| ). A radio point source at 
I = 6.6°, 6 = -0.16° (G 6.6-0 .1) is defined (Altenhoff et 

1983| ), hereafter referred as A83 
A glitching radio pulsar, PSR J1801-23 

s s~-^), lies 
1993|). An up- 



al. |1978i Andrews et al 
in this paper 
(PSR B1758-23, P =416ms, P = 113 x 10 
at the northern radio edge (Kaspi et 



al. 



per limit to this pulsar's characteristic age is estimated at 
58 000 years, and it's spin-down luminosity {E ~ 6.2 x 10'^'* 
erg s~^) is at the lower edge of luminosity values when 
compared to the known gamma-ray (EGRET & COMP- 
TEL) pulsars. 



The age of W28 is estimated (Kaspi et al. 1993) in the 
range 35 000 to 150000 years, with upper and lower limits 
taken from the assumptions that W28 is currently in ei- 
ther the radiative or Sedov phases of expansion. According 



to Kaspi et al. |1993| , the distance of PSR J1801-23 (9 to 
16.5 kpc) derived from it's dispersion measure (DM) ap- 
pears inconsistent with that derived for the remnant. Esti- 
mates for the remnant distance are set at 1.8 kpc (Goudis 



1976 S-D relation) and 3.3 kpc respectively (Lozinskaya 



1981, from mean optical velocities), indicating that the 



pulsar/W28 association is po ssibly a line-of-sight coinci- 
dence. However, Frail et al. (|1993| ) have noted the large 
uncertainty in using the DM as a distance estimate for 
this pulsar due to the high concentration of ionised ma- 
terial in the line of sight, and conclude there is sufficient 
evidence for the pulsar/remnant association. The uniden- 
tified EGRET source 3EG J1800-2338 (95% error cir- 
cle 0.32° radius) (Hartman et al. |l999| ), fisted as 2EG 
J1801-2312 in the second EGRET catalogue (Thompson 



et al. 1996), lies on the edge of the radio shell and was 
thought to be associated with the remnant (Esposito et 
al. |1996i Zhang & Cheng |1998D. 3EG J1800-2338 has a 



relatively hard spectral index (Hartman et al. 1999) with 
no apparent sign of a turnover at 1 GeV (Merck et al. 
1996| ). Lamb & Macomb ( |1997| ) also point out that 3EG 
J1800— 2338 is visible above 1 GeV at 5Aa significance, 
and is centred very close to the A83 radio position. The 
3EG position of the EGRET source is displaced by about 
0.5° relative to the 2EG position, yet still lies comfort- 
ably within the SNR radio shell, and remains a strong 
exam ple of an EGRET source/SNR association (Romero 
et al. 1999 ). The 3EG error circle however, now excludes 
PSR J1801— 23 and the molecular clouds. 



o 
o 
o 

CM 



o 

o 

Q 



' 2EG J1801-2312 




Right Ascension (J2000) 

Fig. 1. Radio continu um ( 327 MHz) for W28, adapted 
from Claussen et al. (1997). Included are the positions 
(and error circles) of the EGRET source 3EG J1800-2238 
(2EGJ1801-2312 in the 2nd EGRET catalogue), OH 
maser sites A to F, PSR J1801— 23 (also the tracking po- 
sition for 1994 data) and th e radio point source, labelled 
A83 (Andrews et al . [l983 ), is the tracking position for 
1995 data. See Sect. \s\ for a discussion. 



W28 lies in a complex region of the galactic plane 
with many HII regions and dense molecular clouds (Woot- 
ten 1981) contributing to the ISM surrounding the SNR. 



Over forty OH (1720 MHz) maser emission sites are con- 
centrated at the eastern and northern edges of the SNR 



(Claussen et al. 1997), along the SNR and molecular cloud 
interface. The distribution of shocked and unshocked gas 
in this region is also consistent with the idea of the SNR 
shock passing through the cloud (Arikawa et al. 1999). 



OH maser emission (1720 MHz) is considered a strong in- 
dicator of coUisional pumping with matter densities ~ 10^ 
cm-3 (Claussen et al. |1997t [19991 ). Enhanced levels of TeV 
7-ray emission via the decay of neutral pions may be ex- 
pected from such areas associate d wit h the masers and 
molecular cloud (Aharonian et al. 1994). Fig. |l| indicates 
the sites of interest in relation to the radio continuum 
emission (327 MHz). The presence of these interesting ob- 
jects make W28 a prime southern hemisphere candidate 
for study at TeV gamma-ray energies. 

We report here on the comprehensive analysis of data 
taken in 1994 and 1995 with the CANGAROO 3.8 metre 
telescope. This work follows analysis of data taken in 1992 



(Kifune et al. 1993) in which weak evidence for a gamma- 



ray signal was reported. At that time, only ON source 
data were collected, making an e stima tion of the back- 
ground rate difficult. Mori et al. ( 1995 ) reported briefly 
on an analysis of 1994 data centred on PSR J1801— 23, in 
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which a ±0.7° field of view was searched. Both ON and 
OFF source data were cohected and no evidence for TeV 
7-ray emission was seen from various point-hke sources in- 
chiding the pulsar and both radio and X-ray maxima. The 
1995 data were centred on the radio position A83, located 
^ 0.3° away from PSR J1801-23. A search for point-hke 
and diffuse sources of TeV emission was carried out on the 
1994 and 1995 datasets out to ±1° from the tracking cen- 
tre of each dataset, using an extended source analysis. We 
have used an improved set of cuts to those used in the anal- 
ysis of data taken on the Vela Pulsar/Nebula (Yoshikoshi 
et al. 1997 ) . These cuts were designed to minimise the loss 
of gamma-ray sensitivity for off-axis sources and in par- 
ticular maintain reliable statistics over the search region. 

2. The CANGAROO 3.8m telescope 

The CANGAROO 3.8 m imaging telescope is located near 
Woomera, Australia (137°47'E, 31°06'S, 160m a.s.L). The 
3.8 metre diameter mirror, of focal length 3.8 metres, is 
used to image the Cerenkov emission from gamma-ray and 
cosmic-ray induced extensive air showers (EAS) onto a 
high resolution multi-phototube camera. The camera con- 
sisted of 224 photomultiplier tubes in 1994. In April 1995 
an extra 32 tubes were added to the corners of the cam- 
era bringing the total to 256 tubes in a 16x16 square grid 
arrangement. However the extra corner tubes have not 
been used in the analysis of the 1995 W28 dataset in or- 
der to retain consistent imaging properties with the 1994 
dataset. Each camera tube is a Hamamatsu R2248 with 
a photocathode size 0.12° x 0.12° on a side, and the to- 
tal field of view of the camera is 2.9°. An event trigger 
is registered when the summed output of triggered tubes 
exceeds a preset threshold, denoted HSUM. Images with 
a minimum of between 3'--^5 tubes, depending on the im- 
age's compactness, trigger the telescope. A vertical event 
rate due to cosmic rays of ~2 Hz is achieved. Monte Carlo 
simulations of the telescope indicate an gamma-ray en- 
ergy threshold of 1.5 TeV at the vertical (Roberts et al. 
1998| ), where the energy threshold is defined as that repre- 



senting the half-maximum of the differential distribution 
of triggered energies. Tracking calibration is performed by 
monitoring the paths of bright (visual mag. 3-6) stars in 
the field of view, providing an absolute tracking accuracy 
of - 0.02° (Yoshikoshi p^96| ). A more detailed technical 
description of the telescope appears in Hara et al. ( 19931 ). 

Data are recorded on clear moonless nights. An ON 
source run is generally followed by an OFF source run 
displaced in right ascension to provide a background run 
of matching zenith and azimuth angle distributions. How- 
ever, since small sections of data from both observation 
seasons were removed due to cloud effects, a normalisa- 
tion, described later, was used in estimating the statistical 
significance of any ON source excess. Pulse charge (ADC) 
and timing (TDC) information for each tube is recorded 
for each event. Calibration of the ADCs and TDCs is 



achieved by recording events triggered with a blue LED 
flasher before each observation. Tube signals are accepted 
as part of an image if they meet a number of criteria: 

1. The TDC value of a tube must lie between ±37. 5ns, 
referenced against the event trigger time. The event 
trigger time is registered when the HSUM threshold is 
met. 

2. The ADC value must be greater than one standard 
deviation above the RMS noise (comprising skynoise 
and electronic noise) for that tube. 

3. The tube must not be isolated. An isolated tube is one 
which is not adjacent to any other accepted tube. 

4. The tube must not have an outlying relative gain value. 
Tubes with relative gains outside the range 1.0±0.3 
contribute significantly to trigger differences across 
the camera face. This factor is particularly important 
when comparing regions over the entire camera face. 

The telescope is an altitude-azimuth type, which intro- 
duces a rotation of the camera relative to the sky about the 
tracking position during data collection. A 'de-rotation' is 
applied to the tube positions to account for this effect, and 
is necessary when considering off- axis sources. The images 
are parameterised according to the moment-based method 
of HiUas ( |1985| ). 

Pre-processing steps designed to minimise the effects 
of electronic interference are described below. The camera 
is divided into groups of eight tubes which share common 
high voltage and other circuitry, and a special cut, box 
(Yoshikoshi |199(: ), is designed to remove images arising 
from electronic contamination, and are concentrated in 
only one or two tube boxes, This box cut, in combination 
with a total ADC sum (adc) cut rejecting events with 
fewer than 200 ADC counts, is very effective at removing 
such artifacts. Monte Carlo simulations show that the box 
cut does not reduce the power of the image cuts, and that 
the optimum adc cut lies at ~200 ADC counts. 

Mirror degradation resulted in an event rate drop by 
about a factor of two from 1994 to 1995, indicating a 
higher energy threshold for the 1995 dataset. The results 
quoted in this work are normalised to a 1.5 TeV thresh- 
old for gamma-rays, using different raw triggering efficien- 
cies for gamma-rays (Sect. ^), which take into account the 
increase in trigger threshold between 1994 and 1995. In 
addition, only events with width >— 0.01 and with the 
number of triggered tubes, ntubes> 4 are accepted. The 
cuts described above are termed noise cuts. 

The ON-OFF statistical significance is calculated fol- 



lowing Li & Ma (1983), before and/or after application of 
all image cuts. 



S = 



ON - P OFF 

y/OWTWOFF 



(1) 



and is used to assess the likelihood of a gamma-ray sig- 
nal. In order to account for the mismatch of observation 
times between ON and OFF source data (and hence zenith 
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angle-dependent event rates), and trigger rate differences 
due to subtle changes in weather conditions and/or tele- 
scope response during observation runs, a normalisation is 
applied to the ON~OFF statistical significance. This nor- 
malisation, P is defined as the ratio of the events available 
for image parametrisation, i.e. after noise cuts. A final sys- 
tematic check on the ON-OFF statistics after application 
of all cuts, performed on a run-by-run basis, is explained 
in Sect. |. 

3. Extended source analysis 

The analysis of extended sources has required the devel- 
opment of new techniques in TeV gamma-ray astronomy 
(Akerlof et al. |l99l|, He ss et al. |1997| , Buckley et al. |1998 , 
Connaugton et al. 199^ ) . These methods use shape param- 
eters such as width and length in combination with source 
position-dependent orientation and location cuts (e.g. al- 
pha, asymmetry and dis). The orientation and location 
parameters are recalculated at every trial source position 
and a skymap is created. In the analysis of CANGAROO 
Vela Pulsar data, Yoshikoshi et al. ( |1997D used a set of 
cuts based on alpha, length, width, concentration and dis. 
A continuous probability distribution was initially used 
to locate the position of the most significant point in the 
skymap. A gamma-ray flux at this point was then esti- 
mated from the ON-OFF excess obtained after using a 
combination of shape and location cuts with alpha< 10°. 
These cuts were optimised using Monte Carlo simulations 
and therefore are a priori decisions. We adopt the same 
a priori philosophy here in order to determine the signifi- 
cance of any result without the need to consider statistical 
penalties. The set of cuts described here were also used in 
the analysis of Centaurus A data described by Rowell et al. 
( 119991 ). The Monte Carlo simulation package, MOCCA92 
(Hillas |1995 ), was used to generate Cerenkov images from 
extensive air showers (EAS). Gamma-ray primaries were 
sampled from a power law above 0.8 TeV with integral 
spectral index —1.6. Cosmic ray primaries, represented by 
a combination of proton, helium and nitrogen primaries, 
were sampled from a power law of spectral index —1.65 
above 1.5 TeV. 

In designing an analysis suited to off-axis and extended 
sources, it is important to consider the off-axis sensitiv- 
ity of the CANGAROO camera, particularly considering 
that it has a relatively small field of view and operates at 
a gamma-ray threshold of about 1.5 TeV. The behaviour 
of various image parameters as a function of gamma-ray 
source position has been investigated. A gamma-ray point 
source was placed at six positions across the camera and 
each position considered independently. The resulting dis- 
tributions of width and miss showed little variation with 
source distance from the camera centre, in contrast to 
those of length and dis (Fig. ||). For increasing source dis- 
tances from the camera centre, larger values of length and 
dis are possible, behaviour which is readily understood 



in terms of camera edge effect reduction on one side as 
the source approaches the camera edge. In particular, the 
length parameter for gamma-rays is very similar to that 
for cosmic-rays for sources beyond about 1.0° from the 
camera centre. Indeed, even for on-axis sources, the field 
of view of the 3.8 m camera imposes some edge effects, 
a consequence of operating at a relatively high gamma- 
ray energy threshold. Thus, for off-axis sources, a reduced 
gamma ray efficiency will result when using cuts derived 
for an on-axis source. 

Firstly, there is a clear case for increasing the gamma- 
ray efRciency while maintaining the quality factor on 
the grounds of improved event statistics. An increase in 
gamma-ray efRciency entails an increase in background ac- 
ceptance, thus improving an estimate of the background 
and diluting any unaccounted systematic effects. It is also 
possible to examine the effect of the gamma-ray efficiency 
on the ON-OFF significance after cut application. We can 
express Eq. (Q) in another way, incorporating the gamma- 
ray efficiency, after application of image cuts: 



2N. 



(2) 



where the gamma-ray signal, (A^^=ON-OFF) and the 
background, Nb {Nb =OFF), represent those prior to im- 
age cuts. The efhciency for gamma-ray selection is given 
by e-y and that for the background (CR), £(,. The quality 
or Q-factor of the cuts is given hy Q = e^/ y/eb . Here, 
we set the normalisation factor from Eq. (|^), f3 — 1. We 
can see that the significance obtained after image cuts is 
dependent on the quality factor and, somewhat slightly, 
on the gamma-ray cut efficiency. For sources with a high 
gamma-ray to CR flux ratio (e.g. >0.1) with Q ~4 (as for 
the CANGAROO telescope), both denominator terms are 
then similar, the sensitivity of S to e-y becomes more ap- 
parent. Such a situation may arise in the case of searches 
for bursts from AGN and/or gamma-ray bursts over short 
time scales. For signal to noise ratios expected of SNRs 
such as W28 however, we would expect only a minor im- 
provement in S from the above arguement. Thus, the main 
motivation for increasing here is simply to work with 
increased statistics. 

The variable cuts on length and for dis were incorpo- 
rated into the cut ensemble. Along with the 3rd moment 
of the image, asymmetry, we use an approximation of the 
distance between the assumed source position and calcu- 
lated source position for each image, D. This distance, 
expressed in units of standard deviation, is given by:- 



D 




dis — dise 



(3) 



where amiss , ^dis are the variances for miss and dis respec- 
tively. The variances, amiss and a^is represent the trans- 
verse and longitudinal errors in the most likely source po- 
sition for an image. D can be characterised as the source 
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Fig. 2. Distributions of length and dis for gamma-rays (solid lines) as a function of gamma-ray source radial distance 
(0.0°, 0.71° and 1.5°) from the camera centre. Cosmic rays (dashed lines) are included for comparison. The co-ordinate 
origin of the camera is shifted to the gamma-ray source position when calculating dis. Noise cuts have been applied. 
In the leftmost plots, real OFF source data, normalised to the simulated CR events are included (dotted error bars), 
as well as mean values for the simulated gamma and CR. 



density function (SDF). When used in combination with 
shape and asymmetry cuts, a cut on D provides a gamma- 
ray acceptance of ~40% for an on-axis point source. D is 
similar to the normalised cluster, or Mahalonobis distance 
for miss and dis (Hillas & West 1991), although here we 



are neglecting cross-term variances. D can also be con- 
sidered a discrete analogue of the probability distribution 
function used by Yoshikoshi et al. ( |1997| ). The longitudinal 
error in source location is greater than the transverse error 
by about a factor of two. Following Akerlof et al. ( 1991] ), 
some reduction in the longitudinal error can be achieved 
by making use of the elongation (defined as length/width) 
of the image such that the expected longitudinal source 
distance is given by disex = 5(1 — (width/ length)) where 
g is an empirically derived constant. For the 3.8 m cam- 
era we derive a value of g = 1.25 using simulations. 
Fig. ^ gives distributions of the parameter D for simu- 
lated gamma-rays and cosmic rays, and a comparison to 



real OFF source data. Linear fits were found for the cut 
on length and value of adis , respectively, as a function of 
source displacement from the camera centre. The total 
cut combination is listed in Table |l|. These cuts (Table |l]) 
provide a constant gamma-ray efficiency of ~40% (and 
cosmic-ray efficiency) and quality factor '^A at the same 
value of D. Without the variable length and dis criteria, 
the gamma-ray acceptance quickly reduces to less than 
30% for sources outside 0.5° from the camera centre for 
no significant improvement in Q-factor. Table ^ gives the 
simulated performance of the cut D for various source po- 
sitions. A comparison to the performance obtained by the 
set o f cuts used in the Vela Pulsar analysis (Yoshikoshi 
et al. 1997 ) is included. These cuts were used to obtain a 
significant excess of gamma-ray-like events from a region 
displaced 0.14° from the pulsar position. The Vela Pulsar 
cuts are based on the same noise cuts listed in Table |], 
except that fixed values of length and dis are used and a 
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Cuts used in this work (see Table 1 ) 
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Table 2. Simulated performance of the image cuts of this 
work (assessed using the quahty factor Q) compared with 
cuts used in the Vela Pulsar analysis. Since the statistical 
errors of Q are ^ 10%, the Vela Pulsar cuts are not statis- 
tically different from those of this work. The Vela Pulsar 
cuts use the same noise cuts as defined in Ta ble W, but 
fixed values of length and dis (Yoshikoshi et al. [l99"7| ). 



Fig. 3. Distributions of D for simulated gamma-rays 
(solid line), simulated cosmic rays (dashed line) and real 
OFF source data (dotted line). 



Parameter Range 



adc > 200 ADC counts 

box < 0.99 

ntubes > 4 

width 0.01<mrf<0.11 

length < 0.32 -|- O.lOd 

asymmetry >0.0 

D < 2.0(7 

adis = 0.21 + Omd 
Cmiss = 0.09 for all d 
dis.. = 1.25(1.0 - ^ 



Table 1. Image cuts used in the extended source analysis. 
See Eq.( ||) for a definition of D. The parameter d is the 
assumed source displacement from the camera centre (in 
degrees). The length cut and value of a dis are dependent 
upon the assumed source position, and are derived from 
linear fits. 



cut on alpha< 10° is substituted for D < 2.0. The Vela 
Pulsar cuts do not provide a constant gamma-ray accep- 
tance over the search region, and the background (cosmic 
ray) acceptance also decreases sharply. 

The resulting flux or upper limit from the search is 
found by dividing the ON-OFF excess, N (or 3a upper 
limit thereof), by the position-dependent exposure (for an 
extended source, averaged over the integration region): 



Fi>E) 



N 



photons cm ^ s ^ 



(4) 



where fj is the raw trigger efficiency for gamma-rays aver- 
aged over the integration region, e-y is a constant gamma- 
ray selection efficiency for the SDF cut, A is the area over 
which gamma-rays are simulated {A — 1.96 x 10^ cm^) 
and T is the total observation time. We set to be the 
average of the simulated gamma-ray cut efficiencies out 
to ±1°, i.e. = 0.41. An unavoidable decline in ry for 
gamma-rays will occur as the source position is displaced 
further off-axis. Simulations show that rj decreases to less 
than half the on-axis value at the corners of the search. 
A linear fit was used to characterise rj at all points within 
±1° such that r]d = ?7o(l-0 — 0.39d) as a function of the 
source displacement from the camera centre, d. We find 
770(1995 data)= 0.12 at the camera centre above the min- 
imum simulated energy of 0.8 TeV for the 1995 dataset. 
The raw gamma-ray trigger efficiency for 1994 data was es- 
timated by scaling the 1995 value by the ratio of observed 
event rates after noise cuts, giving 77o(1994 data)= 0.21. In 
characterising and fj, we are assuming that the gamma- 
ray ffux of an extended source is isotropically distributed. 
To calculate the flux applicable to the energy threshold 
of 1.5 TeV (Roberts et al. 1998) we used the gamma- 
ray spectral index of -1.6 adopted in the simulations. For 
a point-like search, the flux was taken as that from the 
point of interest, using a single value of 77. The nature of 
our gamma-ray selection cuts naturally incorporates the 
gamma-ray point spread function (PSF). The cut on D 
accepts events with derived source positions within an op- 
timal radius. A search for an extended source therefore 
does not require any extra area to account for the PSF in 
addition to the source area itself. In an extended source 
search, at a suitably high number of assumed source posi- 
tions in the region of interest, we sum the events passing 
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1994 Data 



1995 Data 



Time(hrs) 
Images 

After noise cuts 
ON, OFF normalisation (3 
Raw gamma-ray trigger 
efficiency 770 



ON 
20.3 

59909 
37386 



OFF 
23.5 
65483 
40731 



0.92 



0.21 



ON OFF 
37.5 30.0 
97850 76595 
39354 31237 
1.26 

0.12 



all cuts, taking care not to count an event more than once. 
Skymaps of the statistical significance, S, of a gamma- 
ray signal were generated over a ±1° area at 0.05° steps. 
At each grid point, representing an assumed source posi- 
tion, image parameters were calculated and the number of 
events passing the cuts of Table || cumulatively summed 
for all data. Since the resolution of the grid (0.05°), is 
smaller than the effective acceptance area of the cuts (the 
cut on D alone is more powerful than a cut on alpha<l{)°) , 
each skymap point value will not be fully independent of 
its neighbours. 

As a final check on data integrity, the distribution of S 
obtained on a run-by-run basis (i.e. run-by-run skymap) 
was quantitatively assessed for systematic effects. The 
most important systematic effect to consider in this type 
of analysis is the consistency of the trigger threshold over 
the entire camera between ON and OFF source runs. Such 
an effect is difficult to compensate for after the data are 
taken. The Kolmolgorov-Smirnoff (KS) test is used to ex- 
amine the likelihood that the distribution of S of the 
skymap obtained on a run-by-run basis is derived from 
a normal distribution. Over the time scales of a single 
run (<5 hours), we do not expect significant contribu- 
tions from a steady source of TeV gamma-rays of the 
strength expected of a SNR or plerion. A relatively strong 
KS probability of Aa was used to reject pairs of data that 
appeared severely affected by such systematics. For un- 
matched pairs, a well-behaved OFF or ON comparison 
run was used. A total of 10 hours data were rejected using 
the KS test, representing three ON/OFF pairs from the 
1994 dataset. 

4. Results 

Table |^ presents a summary of data accepted for analysis 
from the 1994 and 1995 observing seasons. 

The skymaps of significance S (Fig. ^ for both years' 
data, do not indicate any significant point-like excesses 
over a ±1° search for 7-ray emission. Table ^ summarises 
the 3(7 upper limits for the positions of a number of inter- 



esting sites within the W28 region. It should be pointed 
out that our results will have a systematic error of the 
order ~ 20%, based on uncertainties in the trigger con- 
ditions, mirror reflectivity and spectral index adopted in 



Table 3. Summary of 1994/1995 W28 CANGAROO data 
accepted for analysis. The normalisation factor /3 used in 
adjusting the statistical significance of the post-cut ON- 
OFF excess, and the raw trigger efficiencies for gamma- 
rays above 0.8 TeV (ryo) are included. 



the simulations (Yoshikoshi 1996). For example, in the 
next section we compare our results to a model producing 
a gamma-ray flux with spectral index of ^ —1-1, harder 
than the —1.6 spectrum used in simulations. Such a dif- 
ference however will contribute a systematic of <5%. 
The following positions were considered as potential 



sources; the position A83 given by Andrews et al. (1983), 
the pulsar PSR J1801 — 23, the two strongest (by an or- 
der of magnitude) masers, labelled E & F by Claussen et 
al. ( |1997| ), and the EGRET source 3EG J1800-2338. In 
addition, an extended region of radius 0.25° centred on 
a position to encompass the molecular clouds was con- 
sidered. An average position was used for the masers E 
and F since they are separated by only 0.02°. In assum- 
ing that the EGRET source was point-like, a search for 
the highest point significance within the 95% error circle 
was carried out. Since the statistical degrees of freedom 
for a non-a priori search for point-like emission over the 
skymap is ~100, the highest ON-OFF excess within the 
EGRET error circle must be interpreted with a similar 
statistical penalty in mind. 

The point spread function (PSF) for a pure gamma- 
ray signal can be used to assess the location accuracy of 
the proposed source positions. Monte Carlo simulations 
show that the PSF for gamma-rays increases slightly from 
a standard deviation of 0.2° on-axis to 0.22° for sources 
at the skymap corners. The positions of our candidate 
sources (Table % lie within 1.0° on-axis. Taking a con- 
servative estimate of the PSF as 0.22°, a simplistic esti- 
mate of the source location error is obtained by adding in 
quadrature 0.22/\/l00=0.02° to the the absolute tracking 
precision (0.02°), giving '-^ 0.03°. We have used a value 
of 100 which is representative of the 3cr upper limit ex- 
cesses {N in Eq. (|^)) calculated here. We can therefore 
say that the features listed in Table ^ would be resolved 
by the CANGAROO 3.8m telescope if they were point 
sources of TeV gamma rays. From Table ^ we can also see 
that the positions of the highest significance within the 
EGRET 95% error circle for both year's data are also not 
co-located, with their separation, ^ 0.6°, easily exceeding 
the estimated PSF. 



5. Gamma Ray Production in SNR 

SNR shocks are able to accelerate particles to TeV en- 
ergies. Gamma-rays are produced in secondary reactions 
between these high energy particles and ambient matter 
and radiation fields. The decay of 7r°, produced in ion- 
ion collisions, is the prime hadronic source of gamma-rays 
in SNR. Electrons accelerated to multi-TeV energies in 
SNR give rise to bremsstrahlung and inverse Compton 
(IC) scattering 7-ray production processes in SNR (see 
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W28 1994 



W28 1995 





271 270.5 270 

Right Ascension (J2000 deg) 



269.5 



271 270.5 270 269.5 

Right Ascension (J2000 deg) 



Fig. 4. SDF skymaps of ON-OFF significance for W28 1994 data (left), 1995 data (right). The skymaps extend ±1° 
from the respective tracking centres (indicated by the circled cross) of each year's data. Contour levels indicate the 
same ON-OFF excess over 10 intervals between the maximum and minimum excess values of each skymap. 



1994 DATA 



1995 DATA 



Source 


RA 
(J2000) 


Dec 
(J2000) 


ON 


OFF 


S 


F (>1.5 TeV) 

cm^^s^^ 


ON 


OFF 


S 


F (>1.5 TeV) 

cm^^s^^ 


A83° 


270.21 


-23.34 


644 


698 


+0.1 


<3.36xl0~"^ 


632 


504 


-0.1 


<2. 95x10"'^ 


Mol. Clouds'' 


270.35 


-23.34 


2589 


2843 


-0.3 


<6.64xl0-i2 


2663 


2163 


-0.8 


<5.56xl0-i2 


PSR J1801-23'= 


270.35 


-23.04 


644 


698 


+0.1 


<3. 20x10"^^ 


723 


587 


-0.4 


<3.32xl0~^^ 


Masers E & F"* 


270.47 


-23.31 


765 


822 


+0.3 


<4.14xl0"^^ 


798 


643 


-0.3 


<3.47xl0-i^ 


3EG J1800-2338^ 


270.12 


-23.65 


974 


965 


+2.1 


<8.82xl0"^^ 


804 


575 


+1.9 


<1. 18x10"" 










269.80-^, - 


23.59 




270.43-'^, - 


23.74 


a: Radio point source 


(Andrews ct a 


■ (|l983|) 



















b: Extended source of radius 0.25° encompassing the molecular clouds. 



c: Point source at pulsar position given by from Kaspi et al. (1993). 

d: Point source at average position of 1720 MHz OH mascr sites E and F as defined by Clausscn ct al. (1997). See text, 
e: Highest pointlike significance within 95 % err or circle (0.32° radius) of EGRET source. Quoted position is that for 
the error circle centre (Hartman ct al. 



1999). See also / below. 



f: Positions of each pointlikc maximum significance within the 3EG error circle. 



Table 4. Details of various sites considered as potential emitters of TeV gamma-rays in the W28 region for CANGA- 
ROO 1994 and 1995 data and corresponding 3ct flux upper limits. The statistical significance S is calculated from Eq. 
(|l|) with the normalisation factor P given by Table ^. 



e.g. Mastichiadis & de Jager 1996[ Gaisser et al. 1998D 



The cosmic microwave background (CMB) is usually con- 
sidered the dominant soft photon source with contribu- 
tions (~10%) from the infrared background as seed pho- 
tons for the IC process. The detailed model of Baring et 
al. ( 1999 ) uses these processes collectively (including syn- 



chrotron and bremsstrahlung radiation) to account for ob- 
servations from radio to TeV gamma-ray energies for the 
northern SNR, 1C443. 



Th e mod els of Naito and Takahara (1994) and Drury 



et al. (1994) give predictions of the TeV gamma ray flux 
due to 7r° decay. As a first attempt at explaining the par- 
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tide acceleration processes in W28, we make use of the 
Naito and Takahara model and compare its predictions of 
the 7T° channel to our results. By considering only the 7r° 
production channel, a lower limit on the predicted TeV 
gamma flux can be estimated. A proton parent spectrum 
with differential index —2.1 and exponential cutoff at 100 
TeV has been assumed {dN/dE oc E-"^-^ x exp(-i;/100)), 
i.e. consistent with shock acceleration expected in a SNR. 
The predicted TeV flux will scale according to: 



cx 



Er 



d^ 



(5) 



where Ecr is the energy available for cosmic ray accelera- 
tion, n is the particle number density of the ISM (cm~'^) 
and d is the distance to the SNR (kpc). Ecr is some 
fraction, typically ~10% of the total energy of the SNR 
(canonically lO^-'^erg). In fact, the total SNR energy for 



W28 has been estimated by Rho et al. (1996) at 4x lO^^erg 



from ROSAT X-ray data. 

The most interesting question concerns the possibly 
of enhanced TeV gamma-ray emission from regions of 
high ISM density (Aharonian et al. 1994), where there 
is a greater chance for the interaction of cosmic rays from 
the SNR. The molecular clouds along the northeast and 
northern remnant boundary have been mapped in detail 
by Arikawa et al. ( |1999| ) at the CO J=l-0 and J=3-2 hues. 
The shocked (region undergone passage and compression 
by the SNR shock) component of the clouds is distributed 
along the SNR/clouds boundary, has a mean density of 
> 10'' cm-3, and mass 2 x IO^Mq. The other AxIO^Mq 
of the unshocked gas has a density of < lO^cm"^ and is 
displaced radially outward from the shocked gas by ~1 
arcmin. Clearly, any TeV gamma-ray flux from 7r° decay 
would be dominated by the shocked gas regions, given that 
the unshocked regions of the cloud would have a much 
lower density and lower energy available for CR produc- 
tion, and the mean matter density for regions excluding 
the m olecular cloud is only ^^1.3 cm~'^ (Esposito et al. 
1996 ). Arikawa et al. ( 1999| ) has derived the energy de- 



posited into the shocked gas at 3 x 10 ergs, a value con- 



sistent with the Rho et al. (1998) estimate for the total 
SNR energy when considering the volume filling factor V, 
between the clouds and the SNR. V is simplistically esti- 
mated at < 0.01, by taking the mass and density of the 
shocked gas (given above) , assuming the cloud consists of 
H2, and a value of lOpc for the SNR radius. We also as- 
sume that ~10% of the available SNR energy goes into 
cosmic ray production, a reasonable value for a Sedov- 
phase SNR, and also consistent with the measured ener- 
getics for the SNR and clouds. Thus we adopt values of 
3x 10^'^ergs for Ecr, and 10** cm^'^ for n in Eq. (||). A work- 
ing band on our flux prediction is obtained if we assume 
a range of values for d from 1.8 to 3.3 kpc, as discussed 
m Sect. I. In Fig. |, we compare the model predictions 
based on the above scalings to our upper limit obtained 
from an extended source of radius 0.25° encompassing the 



10" 



10 



-6 



o 10" 
o 



10 

M 

A 



-10 



10' 

10" 



-12 



-14 
10" 




HEG .J1801-2312 (Merck etal. 1996) 



E„=3xl0*' ergs'-;^, 

11= 1x10" cm 

d=L8 to 3.3 kpc ' 



^0 



.2 



10" 10^ 10 

7-ray Energy (GeV) 



4 



10"- 



Fig. 5. Comparison with our upper limit and the EGRET 
flux of 3EG J1800-2338 (Merck et al. 1996) with a 
model predicting the TeV gamma-ray flux due to tt° de- 



cay (hashed area, Naito and Takahara 1994). The hashed 
region is bounded by limits on the predicted flux when 
assuming published ranges of values for various scaling 
parameters (eq. ||) deflned for the regio n of shocked gas 



19991) . This region is 



surrounding W28 (Arikawa et al. 
expected to dominate those that would produce a TeV 
gamma-ray flux from the decay of tt°. Our upper is for 
an extended source containing the molecular clouds (see 
text). The model is practically upper-bounded by the flux 
from 2EG J1801— 2312, and parent-proton energies are de- 
rived from a power- law of index —2.1 (differential) with 
an exponential cutoff at lO^^'cV. 



clouds from 1994 data (the highest of our upper limits). 
The flux from 2EG J1801-2312 (essentially the same as 
3EG J1800-2338) and its straight extrapolation to TeV 
energies is also included as any 7r° gamma-ray flux will be 
limited by the EGRET measurement. 

Our upper limit lies an order of magnitude below the 
straight extrapolation of the flux from EGRET (dashed 
line of Fig ||) , and is able to place some constraint on the 
prediction of a n° gamma-ray flux from the shocked gas 
region. In order to accommodate our upper limit however 
a cutoff and/or slighter steeper parent spectrum than -2.1 
appears necessary. In assuming an exponential cutoff at 
100 TeV for the parent spectrum of accelerated hadrons, 
we assume that W28 follows the 'standard' picture of par- 
ticle acceleration in SNR, with particle energies limited 
by radiative losses, flnite age of the shock and particle es- 
cape. When the cutoff for electrons is due only to radiative 
losses however, it can be expected that the hadron spec- 
trum will continue (Reynolds & Keohane 199E ). Apart 
from the long-established 'knee' at '--^ 5 x 10^^ eV in the 
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all-particle cosniic ray spectrum there is now direct ex- 
perimental evidence pointing to a continuation of the pro- 
ton/helium spectra up to at least ~ 800 TcV (Asaki- 
mori et al. 199S), implying that strong cutoffs may not 
be required. Our result here is consistent however, with 
previous comparisons of upper limits (from other SNRs) 
to hadron-induced gamma-ray models (e.g. Buckley et al. 
1998| , Allen et al. |1995| and Prosch et al. |1996D which seem 



cations of our results in relation to the origin of galactic 
cosmic rays. 

Further data on W28 will no doubt be taken with the 
recently completed CANGAROO II telescope (Yoshikoshi 
et al. 1999), which, at the very least will provide tighter 



to require some cutoff below the knee energy/and or spec- 
tra steeper than -2.1. Further constraints may arise if elec- 
tronic components are considered. Particularly, electronic 
bremsstrahlung may dominate over the inverse Compton 
component due to the very high density of target matter 
in the clouds. 

The above discussion aside, the location of 3EG 
J1800— 2338 by itself, makes it's interpretation difficult in 
terms of simple CR-matter interaction, or as the re sult o f 
a pulsar-powered process (proposed by Merck et al. 1996| ) . 
Many of the promising sites within W28 for gamma-ray 
production are now outside the 95% error circle, leaving 
just the filled X-ray centre, and southern/ western portions 
of the radio shell. At best, we are in a position to rule out 
the interpretation of the EGRET source as resulting to- 
tally from 7r° decay gamma-rays with an unlimited parent 
spectrum, and to place limits on the parent spectral in- 
dex/cutoff energy combination. 

6. Conclusion 

A search for TeV gamma-ray emission from the south- 
ern SNR W28 was carried out by the CANGAROO over 
two observation seasons (1994 and 1995) using the atmo- 
spheric Cerenkov imaging technique. An analysis provid- 
ing a consistent gamma-ray acceptance and quality factor 
for extended sources was used. A number of sites within 
a search region of ±1° were considered as potential point- 
like and diffuse sources of TeV gamma-ray emission. No 
evidence was found for the emission of TeV gamma-rays at 
any of these sites which include those from the strongest 
two masers, an EGRET source, a radio pulsar (all as point 
sources) and a diffuse region containing the molecular 
clouds that appear to be interacting with the remnant. 
Our 3(7 upper limit from this diffuse region at 6.64xl0~^^ 
cm-^s-i > 1.5 TeV, and the flux of the EGRET source 
3EG J1800— 2338 were compared with gamma-ray flux 
predictions from the model of Naito and Takahara (1994). 
Under this framework, our upper limit rules out a straight 
extrapolation of the EGRET flux to TeV energies. It also 
constrains somewhat the flux expected from the shocked 
region of gas in the molecular cloud, placing limits on the 
parent spectra for hadrons and/or cutoff energy. Our re- 
sults suggest that the EGRET source probably does not 
result entirely of tt° gamma-rays. This fact is supported by 
it's location in relation to the molecular clouds. In a later 
paper, we will consider electronic bremsstrahlung and in- 
verse Compton scattering and discuss the broader impli- 



constraints on models of gamma-ray production for this 
interesting source. 
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